Bivariate hierarchical Hermite spline quasi--interpolation

Spline quasi-interpolation (QI) is a general and powerful approach for the construction of low cost and accurate approximations of a given function. In order to provide an efficient adaptive approximation scheme in the bivariate setting, we consider quasi-interpolation in hierarchical spline spaces. In particular, we study and experiment the features of the hierarchical extension of the tensor-product formulation of the Hermite BS quasi-interpolation scheme. The convergence properties of this hierarchical operator, suitably defined in terms of truncated hierarchical B-spline bases, are analyzed. A selection of numerical examples is presented to compare the performances of the hierarchical and tensor-product versions of the scheme

Discontinuity Detection by Null Rules for Adaptive Surface Reconstruction

We present a discontinuity detection method based on the so-called null rules, computed as a vector in the null space of certain collocation matrices. These rules are used as weights in a linear combination of function evaluations to indicate the local behavior of the function itself. By analyzing the asymptotic properties of the rules, we introduce two indicators (one for discontinuities of the function and one for discontinuities of its gradient) by locally computing just one rule. This leads to an efficient and reliable scheme, which allows us to effectively detect and classify points close to discontinuities. We then show how this information can be suitably combined with adaptive approximation methods based on hierarchical spline spaces in the reconstruction process of surfaces with discontinuities. The considered adaptive methods exploit the ability of the hierarchical spaces to be locally refined, and fault detection is a natural way to guide the refinement with low computational cost. A selection of test cases is presented to show the effectiveness of our approach